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Slater and Dunkl reported formulas-involving (generalized) hypergeometric functions 
(p-Fp_i)-for the Hilbert-Schmidt moments of two sets of univariate probability distributions 
pertaining to the entanglement/separability of 2 x 2 quantum systems [1, p. 30] (cf. [2]). 
One ("balanced") set of (determinantal) moments had the form (p denoting a 4 x 4 density 
matrix, and p'^^, its partial transpose) 

^ (2n)!(l + a),Jl + 2a),„ ^ / -n, a, a + i, -4n - 1 - 5a _ \ 

and the other set of ("unbalanced") determinantal moments, the form 

nPT\n^ ^ n! (g + 1) J2a + 1)^ 

, (-2n-l-5a)Ja)„(a+|)^ ^ / -i^, -n, « + 1, 2« + 1 _\ 

24^^ (3a+ 1)^ (6a+ ^ ^ U - ^, « + 2 + 5a, 1 - n - a, | - n - a' / 

(Because of the denominator parameter 1 — n it is necessary to replace the 5F4-sum by 1 
to obtain the correct value when n = 1. Both hypergeometric functions are terminating (as 
well as balanced in the Pfaff-Saalschutzian sense [31 sec. 2.2] in nature.) Here a is a Dyson- 
index- like parameter that takes the value |, 1 and 2 for real, complex and quaternionic 
systems, respectively. 

The probability distributions having the first set of balanced moments extend over the 
range [— 2~^^ ■ 3~^,2~^^], with the nonnegative interval [0,2^^^] corresponding to separable 
systems. The probability distributions possessing the second set of unbalanced moments 
extend over [— 2~^,2~®], with the nonnegative interval [0, 2~®] corresponding to separable 
systems. 

We are interested in finding the exact probability distributions yielding these sets of 
moments-the accomplishment of which task would presumably yield further insight into the 
nature and character of the derived separability probabilities. In undertaking such an effort, 
to begin, we follow-as best we can-the analytical framework, involving Mellin transforms 
and Meijer G functions, employed by Penson and Zyczkowski in their study of Fuss-Catalan 
and Raney distributions [1]. 

We have largely worked with the second set of (unbalanced) moments in recent efforts 
of ours [51 Ej. In fact, in [^, making use of this particular set, we were able to report the 
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following "concise" formula for the probability -P(«) that a 2 x 2 system is separable (that 
is, the cumulative probability over [either of] the indicated nonnegative intervals) 



P(«) = E-o/(« + 0, (1) 

where 

ff ^ P( ^ P( g(a)2-^"-6r(3a+|)r(5a + 2) 

f{a) = P{a)-P{a + l)= ^'itrf^ ^ i3V 

3r(a + l)T{2a + 3)r(5a + y) 



and 



q{a) = 185000a^ + 779750a^ + 1289125a^ + 1042015a2 + 410694a + 63000 = (3) 

a(5a(25a(2a(740a + 3119) + 10313) + 208403) + 410694) + 63000. 

It was, then, concluded (using strongly convincing numerical evidence consisting of expan- 
sions of hundreds of digits) that for the (nine-dimensional) two-rebit systems, P(| 
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for the (fifteen-dimensional) two-qubit systems, P(l) = and for the (twenty-seven- 
dimensional) two-quater(nionic)-bit systems, P(2) = 

Nevertheless, despite these substantial successes with the particular use of the second 
set (convergence of probability-distribution reconstruction procedures p] being much slower 
with the first set), it appears that the first set of moments is more immediately amenable to 
the implementation of the Penson-Zyzckowski approach for reconstruction of the complete 
probability distributions of interest. 

Since (quite significantly, it would appear) the ranges of the indicated probability distri- 
butions include negative-valued intervals, it appears necessary-for application of the Mellin 
transform-to appropriately modify the distributions to eliminate nonnegative domains. If 
we do linearly transform the balanced moments of the probability distributions defined over 
[— 2~^^ ■3~'^, 2^^^] to the moments of probability distributions defined over the (nonnegative) 
unit interval [0,1], we obtain for the new n-th moments (using the binomial expansion) 

m ip""ir)io,i - ( 2-1. . 3-3 + 2-., )>^^. Q ■ 3')- m Ip""!)-") ■ w 

(Of course, it would be highly desirable to have a more explicit expression for 
((IpI |p^'^|)")[o 1] than this one-a task we are pursuing.) Now, (the original, untrans- 
formed) ((|p| Ip-^-^l)""™) moment occurring in the right-hand side can be expressed-using 
the expansion formula for hypergeometric functions, as well as the Gauss multiplication 
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(2 3 le { m+n) 12a Gamma \ — + m - n] Gamma [Jt + m - n] Gamma \ m + n] Gamma [1 - m + n] 

r 1 k 5 a-, r k 5 a-, r 1 k 5 a-, 

Gamma + — + m - n Gamma — + m - n - Gamma — +— +m-n - 

L44 4 i I 4 4 i ^4 4 4 ^ 



Gamma [m - n - a] Gamma ^ — + m - n - a| Gamma ^m - n - — j 



r 1 k 5 a- 

Gamma + +m-n - 

L 2 4 4 

Gamma — + m - n - — Gamma m + n + — Gamma \ 1 - m + n + — Gamma [ Jc + a ] 

L 2 2^ ^2 2^ ^ 2^ 

Gamma — + Jc + a Gamma m + n + a Gamma [1-m + n + a] Gamma — + 3 a Gamma \ — -t- 6 a\ 

I 2 J ^ 2 2 2 1' 

I r 1 -| rl Sa-, r SO-, 

Gamma [1 + k] Gamma [m - n] Gamma — + Jc + m - n Gamma + m - n - Gamma m - n - 

\ 2 ^ ^ 4 4 ^ ^ 4 ^ 

rl Sa-, r 1 5 a ^ r k 1 r 1 k 1 

Gamma — + m - n - Gamma — + m - n - Gamma — + m - n - a Gamma — + — + m- n- a 



L2 4 i 2 J 2 2 J 

rJC a-, rlJC a-, rl-. 

Gamma — + m - n - — Gamma — + — + m- n- — Gamma [ a ] Gamma — + a 

y 2 2^ ^22 2^ ^ 2 ^ 

r 5 3 a-, V 3 3 a-, r7 3 a^ 

Gamma [1 + a] Gamma m + n + Gamma m + n + Gamma m + n + 

Lg ^4 2-' 2 ^ 

r5 3a-, r5 3a-| r 11 3 a-, \ 

Gamma m + n + Gamma m + n + Gamma m + n + Gamma [1 + 2 a] 

^8 2 ^ ^4 2 ^ ^8 2 ^ ] 

FIG. 1: Quantity which when summed over k from to n — m yields Ip-^^])""™^ 

theorem (for subsequent use of the inverse Melhn transform)-as the sum over k from 
to n — m of the quantity in Fig. [l} (Fig. [2] gives the function-setting a = 1-which 
when summed over k from to n, then added to (2^^ ■ 3'^)"", with the resultant 
sum multiphed by ( 2-12.3-3+2-16 yields the unit-interval moments ((|p| |p'^^|)")[o j^]-) 
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1 r 3 T 
2^-" ° Gamma I - + k I 

Dlf f erenceRoot ^Function ^ |y, n|, ||-6 + 4n + k- 4nj |-5 + 4n + k- 4nj |-4 + 4n + k- 4nj 

(-3 + 4n + k-4n) (n + k - n) (7-4n + 4n) (9 - 8 n + 8 n) (ll - 8 n + 8 n) (l3 - 8 n + 8 n) 
(15 - 8 n + 8 n) Y[n] + (-48 648 600 n - 2 883 584 n" - 48 648 600 k + 46 216 170 k^ - 

16 081 065 k^ + 2 432 430 k* - 135 135 k^ - 65 536 (379 + 196 k - 332 n) + 48 648 600 n - 

392 889 906 k n + 298 672 137 k^ n - 86 765 970 k^ n + 12 154 455 k^ n - 948 516 k^ n + 

346 176 072 - 1 368 013 392 k + 838 396 260 k^ - 207 992 824 k^ + 28 098 660 k' - 

2 089 088 k= + 1 082 685 168 - 2 720 578 752 k + 1 355 060 736 k^ - 

292 540 096 + 35 634 304 k* - 1828 352 k= + 1960 350 528 - 3 443 835 264 k n* + 

1400 219 520 k^ - 256 893 952 k^ n* + 23 414 272 k* n* - 569 344 k= + 2 277 476 608 n= - 

2 911 979 520 k n° + 948 948 992 k^ n° - 132 038 656 k^ n° + 6 434 816 k^ - 16 384 k^ + 

1772 422 144 - 1 654 177 792 k + 399 245 312 k^ - 32 210 944 k^ + 

278 528 k* n' + 926 973 952 n' - 600 113 152 k n' + 87 162 880 k^ n' - 1 835 008 k^ n' + 

311 099 392 n° - 119 275 520 k + 5 767 168 k^ + 58 458 112 - 8 388 608 k n' + 

4 194 304 - 65 536 (-3829 + 335 k + 290 k^ - 36 (109 + 40 k) n + 1008 n^) - 

32 768 n' (26 858 + k (-13 103 + k (-607 + 376 k) ) + 63 100 n - 

2 k (4165 + 1942 k) n + 72 (493 + 124 k) - 2784 n^) + 
64 n* (30 723 789 + 2 k (-26 900 643 + k (10 559 919 + 4 k (-392 234 +k (9659 + 1696 k) ) ) ) + 

178 543 028 n + 32 k (-7 012 557 + 2 k (1025 653 + k (-79 532 + k (-6293 + 412 k)))) n- 

240 (-1 725 409 + 16 k (95 239 + k (-17 222 + k (-338 + 145 k) ) ) ) + 

2560 (192 299 + 2 k (-53 677 + 4 k (694 + 167 k) ) ) - 

35 840 (-8576 + k (2155 + 202 k) ) n* + 1032 192 (83 + 2 k) + 4 472 832 n') + 
4 (86 544 018 + k (-342 127 764 + k (209 661273 + k (-51 624 958 + 

(6 589 209 - 397 856 k) k) ) ) + 812 511 540 n - 48 k (42 545 331 + 
k (-21113 208 + k (4 399 867 + 2 k (-227 849 + 6508 k) ) ) ) n + 96 (30 677 133 + 
2 k (-26 902 803 + k (10 749 567 + 4 k (- 446 990 + k (27 695 + 292 k) ) ) ) ) + 
128 (44 558 861 + 8 k (-7 060 941 + 2 k (1092 019 + k (-120 356 + k (781 + 196 k) ) ) ) ) 

- 3840 (-1728 145 + 16 k (98 101 + k (-20 795 + 814 k + 64 k^) ) ) + 
24 576 (195161 + 2k (-58 888 + k (6043 + 236 k) ) ) n° - 

114 688 (-18 070 + k (5795 + 26 k) ) n' - 1179 648 (-389 + 28 k) n' + 31260 672 n') - 
n (-3 k (130 963 302 + k (-99 557 379 + k (28 921 990 + k (-4 051 485 + 316 172 k) ) ) ) - 
8 (-86 544 018 + k (342 065 556 + k (-209 630169 + k (51811 582 + 

k (-6 806 937 + 460 064 k) ) ) ) ) n - 48 (-67 688 559 + 4 k (42 527 187 + 
k (-21 143 016 + k (4 485 403 - 506 242 k + 20 792 k^) ) ) ) - 256 (-30 653 805 + 
2 k (26 903 883 + k (-10 844 391 + 4 k (474 368 + k (- 36 713 + 410 k) ) ) ) ) + 
256 (44 520 413 + 8 k (-7 085 133 + 4 k (562 601 + k (-70 384 + k (2159 + 44 k) ) ) ) ) n* - 
6144 (-1729 513 + 8 k (199 064 + k (-45 163 +k (2780 + 47 k) ) ) ) n= + 32 768 (196 592 + 

k (-122 987 +k (15 353 + 40 k) ) ) + 65 536 (37 058 + k (- 13 075 + 326 k) ) n' - 
589 824 (-835 + 92 k) n" + 34 865 152 n') - 1024 (-1 722 673 - 6 061984 n + 
16 (k (92 377 + k (-13 649 + 2 k (-7 45 + 113 k) ) ) + 4 k (48 466 + (491 - 1100 k) k) n + 
28 (-16 234 + k (2825 + 782 k) ) - 7392 (25 + 4 k) + 672 n*)) - 256 n= 
(8 934 841 + 41376 984 n + 8 (k (-1397 673 + 4 k (99 247 +k (-5912 + k (-983 + 52 k) ) ) ) + 
24 k (-187 616 + k (30871 + (1828 - 371 k) k) ) n + 48 (190868 + 

k (-102 143 + k (2285 + 1768 k) ) ) - 224 (-33 386 + k (7135 + 1186 k) ) + 
4032 (607 + 52 k) n* + 80 640n=)) - 16 (67 730 031 + 491207 376 n + 
4 (k (-42 563 475 + k (21 083 400 + k (-4 314 331 + 405 154 k - 5240 k^) ) ) + 

8 k (-26 901723 + k (10 654 7 43 + 4 k (-419 612 + k (18 677 + 994 k) ) ) ) n + 

8 (44 597 309 + 8k (-7 036749 + 8k (264 7 09 + k (-24 986 + k (-689 + 76 k) )) )) - 

320 (-1726 777 + 8 k (193 340 + k (-38 017 +k (476 + 209 k) ) ) ) + 

2560 (193 730 + k (-112 565 + k (8819 + 904 k) ) ) - 
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7168 (-35 222 + 215 k (47 + 2 k) ) n= - 86 016 (-721 + 20 k) + 3 883 008 n') ) ) 
y[l + n]+27 648 (l+n) (-2+2n + k-2n) (- 1 + 2 n + k - 2 n) ^ (2 n + k - 2 n) 
(l + 2n + 2k-2n) (-n + n) (l-n + n) 
(3-2n + 2n) (3-4n + 4n)y[2+n] ==0, 
y [0] ~ 0, y [1] ~ (Gamma [2 (1 + n) ] Gamma [1 + 2 n] Pochhammer [3, 2 n] 

Pochhammer [ - 6 - 4 n , k] Pochhammer [-n, k])y^ 
|pochhanimer ^ — , 2 nj Pochhammer ^ — r 4 n j Pochhammer [-l-2n, k] 

Pochhammer ^ — - n, kj Pochhammer [- 2 (1 + n) , |] ] ' 
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FIG. 2: The function-solving the indicated hnear difference equation-which when summed over 

^2-12.3-3+2- 



k from to n, then added to (2^^ • 3^) with the resultant sum multiplied by ( 0-12.0-3 , 0-16 )"'; 



yields the transformed (unit-interval) moments ((|p| |p^"^|)")[o 1] two-qubit case a = 1. 

I , r- I ^ 1 r^i 

I 2 3 16 { 1 m) 12a ^ ^ Gamma[k + a] Gamma — + Jfc + a Gamma — + 3 a Gamma — + 5a 

[ L 2 J L 2 J L 2 J 

.1 1 1 1 

MeijerG\\\ m, - k - m, — {-l-2m+a), — {-2m + a), m + a, -m + a, 

III 2 2 2 2 

1 1 11 



{- 2 - k - 4 m + 5 a) , — {- 1 - k - 4 m + 5 a) , — {-k-4m+5a), — {1 - k - 4 m + 5 a) \ , 

A A ^ J 

^ ^ 11 

, — { - 1 - 4 m + 6 a) , — {l-4m+6a)\\, 
4 A J J 

m, -m, — {-l-2m+a), — {-2m + a), m + a, -m + a\, 

o o o > 



4 4 4 4 

3 3a 1 3a 1 3a 3 3a 1 1 

jn + , m + , — - m + , — - jn 

8 28 28 28 24 4 

11 1 1 



2 

1 1 1 1 k k 

■ k - m, -m, — {-l-k-2m+a), — {-k-2m+a), m+a, m+a, 

2 2 2 2 2 2 

5a 1 



-m + — 
4 



-, - {-2-4m+5a), - {-l-4m+5a), - {1 - 4m+ 5a)]], 
4 4 4 ' ' 



1 Gamma \ 1 + k] Gamma fal Gamma | — + a| Gamma [1 + a] Gamma [1 + 2 a] 



FIG. 3: Inverse Mellin transform of quantity in Fig. [T] 

Applying the inverse Mellin transform to the quantity in Fig. [ij having made the 
required transformation n cr — 1, we obtain the expression in Fig. [3j Then, 
Fig. |4] shows an equivalent form-in terms of hypergeometric functions-to the Meijer 
G term in Fig. [s] with the summation index m having been set to zero (cf. [1]). 
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FIG. 4: Equivalent form-in terms of hypergeometric formulas-to the Meijer G function term in 
Fig. [3] with the summation index m having been set to zero. 

So, certainly highly formidable obstacles remain to the achieving of our stated goal of 
explicitly constructing the exact univariate probability distributions-parameterized by the 
Dyson-index-like parameter a-for 2x2 quantum systems that yield the (balanced) univariate 
Hilbert- Schmidt determinantal moments ((|p| |p'^"^|)")obtained by Slater and Dunkl [1]. 

One immediate issue to be addressed-so that Mellin transform methods might be more 
readily employed-is the development of explicit formulas for the determinantal moments, 
transformed-to avoid negative domains-so that they correspond to probability distributions 
over the range [0,1]. Perhaps it is possible to exploit the fact that the hypergeometric 
functions in both sets of moments presented at the outset of the paper, are balanced (in the 
Pfaff-Saalschutzian sense [31 sec. 2.2] ) and terminating in character. Further, the use of 
Euler's integral representation of the generalized hypergeometric function [2| might prove of 
value. 

As another approach to ascertaining the properties of the probability distribution func- 
tions in question, we have used the Mathematica-implemented Legendre-polynomial-based 
reconstruction algorithm of Provost |E], that we have previously applied with considerable 
success (HE], to determining the intercepts of the Hilbert-Schmidt determinantal proba- 
bility density functions. That is, we seek the values of these probability density functions 
at which (the x-variable) \p^'^\ is zero. In Fig. [s] we show the seventy y-intercepts, as a 
function of the Dyson-index-like parameter a = \ (rebits), 1 (qubits), |,2 (quaterbits), 
. . . , 35, for the class of probability distributions, extending over [— 2~^^ ■ 3^^, 2~^^] based on 
the first (balanced) set of moments. (The first 1900 determinantal moments were employed.) 
In Fig. |6] we show the seventy y-intercepts, as a function of the Dyson-index-like parame- 
ter a = ^, 1, |, . . . , 35, for the class of probability distributions, extending over [—2~^, 2~^] 
based on the second (unbalanced) set of moments. (The first 2000 determinantal moments 
were employed.) We hope to be able to discern (increasing further the numbers of moments 
employed) exact values for these y-intercepts, which would then hopefully cast light on the 
specific nature of the Hilbert-Schmidt determinantal probability density functions that have 
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y-intercept 
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1.0x10^ 
8.0x10^ 



■•••-^-'lO 15 20 25 30 35 



a 



FIG. 5: The y-intercepts-at which Ip'^^l = 0-as a function of the seventy values of the Dyson- 
index-hke parameter a = ^,1,|,...,35, for the class of probability distributions, extending over 
[— 2~^^ • 3~^,2~^^] based on the first (balanced) set of moments. The first 1900 determinantal 
moments were employed in the probability-distribution Legendre-polynomial-based reconstruction 
process 



been the subject of this communication. 

Another possible use of the Legendre-polynomial-based probability-density reconstruction 
process |8j might be to determine the modes of the yet-unknown separability probability 
density functions as functions of a. Along such lines, in Fig. [7| we present-based on the first 
1250 unbalanced moments, setting a = 1, an approximation to the two-qubit probability 
distribution defined over [—j^, (An estimate-based upon the first 500 moments-of 

the median of the distribution is \p^'^\ = —0.00691863.) The cumulative (separability) 
probability over the nonnegative interval [0, appears to be equal to J^, as we have 
recently argued PQ El E] (and indicated at the outset of this paper). In Fig. [sj we show the 
two-rebit counterpart [a = ^) to this curve, and in Fig. |9] the two-quaterbit counterpart 
{a = 2). (Estimates of the medians-based upon the first 500 moments-of the last two 
distributions are, respectively, \p^'^\ = —0.00562687 and |p^^| = —0.0121435.) 

For the two-rebit (a = ^) systems, we have the (univariate Hilbert-Schmidt determinantal 
moment) formulas [9][eq. (3.2)] eq. (1)] (cf. [lOl Theorem 4]): 

/I n..^3-2fcr(2A; + 2)r(2A; + 4)\ . , 

m=H^''' r(4feVio) )■ 

Now, IpI G [0, It might be analytically convenient (for convolution purposes,. . . [cf. [TTl 
App. B]]), at some point, to linearly transform this variable to possess the same range as 
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FIG. 6: The y-intercepts-at which \p^'^\ = 0-as a function of the seventy values of the Dyson- 
index-like parameter a = 35, for the class of probability distributions, extending over 
[— 2^^,2~®] based on the second (unbalanced) set of moments. The first 2000 determinantal mo- 
ments were employed in the probability-distribution Legendre-polynomial-based reconstruction 
process 
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FIG. 7: Approximation-based on the first 1250 unbalanced moments, having set a = 1-to the 
two-qubit separability probability distribution over 



.1. J_] 

16' 256 J' 



\p^'^\, that is, [—jq, ^] (as opposed to transforming the variable |p^^| to range over the 
unit interval). Under such a transformation, the new forms of the moments (j5]) are given 



by the function in Fig. [10} solving the indicated linear difference equation. Further, the 
probability distribution function (/(y), with y = 17\p\ — ^) yielding these moments is (cf. 
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FIG. 8: Approximation-based on the first 1250 unbalanced moments, having set a = g-to the 
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FIG. 9: Approximation-based on the first 1250 unbalanced moments, having set a = 2-to the 
two-quaterbit separability probability distribution over [— j^, jm]- 
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Dif ferenceRoot |^ 

Function [{y, n} , |-n (l + n) (2 + n) y[n] +8 (l + n) (2 + n) (l29 + 26n) y[l + n] + (362 496 + 368 640 n h 

116 736 li^ + 11520 h^) y[2 + h] + (13 234 176 + 9 902 080 h + 2 414 592 + 192 512 n^) y[3 +n] + 

(130 744 320 + 78 577 664 h + 15 728 640 + 1048 576 h^) y[4 + n] =-- 0, 

63 453 39227 

y[l] == 1, y[2] , y[3] == , y[4] } [k] 

1144 146 432 222 576 640 ■' 



FIG. 10: Formula satisfied by the two-rebit (a = |) Hilbert-Schmidt moments after their trans- 
formation so that the original range [0, ^gg] of -|:^e determinant of the density matrix matches the 
range [—jq, ^] of the determinant of the partial transpose. 



